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We consider a wide class of static spherically symmetric black holes of arbitrary dimension with a 
photon sphere (a hypersurface on which a massless particle can orbit the black hole on unstable circu- 
lar null geodesies). This class includes various spacetimes of physical interest such as Schwarzschild, 
Schwarzschild-Tangherlini and Reissner-Nordstrom black holes, the canonical acoustic black hole or 
the Schwarzschild-de Sitter black hole. For this class of black holes, we provide general analytical 
expressions for the Regge poles of the S-matrix associated with a massless scalar field theory. This 
is achieved by using third-order WKB approximations to solve the associated radial wave equation. 
These results permit us to obtain analytically the nonlinear dispersion relation and the damping 
of the "surface waves" lying close to the photon sphere as well as, from Bohr-Sommerfeld-type 
resonance conditions, formulas beyond the leading-order terms for the complex frequencies corre- 
sponding to the weakly damped quasinormal modes. 
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I. INTRODUCTION 

Quasinormal modes (QNMs) of black holes (BHs) have 
been studied for nearly 40 years due to their importance 
in the context of gravitational wave astronomy. In the 
last decade, there has been moreover an increase of activ- 
ity in BH QNM studies motivated by potential applica- 
tions in analog models of gravity, quantum gravity, string 
theory and related topics (TeV-scale gravity, AdS/CFT 
correspondence, alternative theories of gravity, BH area 
quantization, phase transitions in BH systems,...). For 
excellent reviews on the status of QNMs prior to 1999 
and on their relevance to gravitational wave astronomy, 
we refer to the articles by Kokkotas and Schmidt |l| and 
by Nollert 0. For a more recent review on BH QNMs, we 
refer to the article by Berti, Cardoso and Starinets Q: 
it updates the two previously cited articles and it also 
presents the aspects of QNM physics linked to gauge- 
gravity duality; it includes, furthermore, an interesting 
historical introduction on the subject as well as a useful 
impressive bibliography on all the aspects of BH physics 
linked to QNMs. 

Immediately after the publication of one of the first 
papers on QNMs by Press [1] where he identified the 
gravitational ringing of the Schwarzschild BH as due 
to its "free oscillations", Goebel suggested a physi- 
cally intuitive interpretation of the associated QNMs 
they could be interpreted in terms of gravitational 
waves in spiral orbits close to the unstable circular pho- 
ton/graviton orbit at r = 3M which decay by radiating 
away energy (here M denotes the mass of the BH). This 
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appealing interpretation has been developed by other 
authors for various field theories denned on BH back- 
grounds using the eikonal approximation, i.e., in a frame- 
work based on geodesies and bundle of geometrical rays 
(see Refs. 0-GJ] as well as Ref. [14[ for a more mathemati- 
cal approach) . It has permitted them to obtain analytical 
approximations for the leading-order terms of the charac- 
teristic complex frequencies of various BH spectra from 
an interpretation in terms of massless particles "trapped" 
near unstable circular null geodesies (see, more particu- 
larly, Ref. [ll[ where the relation with the Lyapunov ex- 
ponent corresponding to geodesic motion is clearly em- 
phasized). 

A potentially much richer implementation of the 
Goebel interpretation of BH QNMs which is not lim- 
ited to purely geometrical considerations but based on 
wave/field theory and which goes beyond the leading 
oder terms has also been formulated (l5l - [r7j (see also 
Ref. [l2j ). It uses complex angular momentum (CAM) 
techniques (or, in other words, the Regge pole machin- 
ery) which play a central role in scattering theory. Since, 
as noted by Chandrasekhar and coworkers [l8| (see also 
Ref. [13), BH perturbation theory can be formulated as a 
resonant scattering problem, CAM techniques arise natu- 
rally in BH physics. For reviews of the CAM method, we 
refer to the monographs of Newton [2(| , Nussenzveig [2lT ] 
and Collins [22j as well as to references therein for vari- 
ous applications in quantum mechanics, nuclear physics, 
high energy physics, electromagnetism and seismology. 

Some years ago, the CAM method was used in grav- 
itational wave physics by Chandrasekar and Ferrari (23j 
to express the flow of energy due to nonradial oscillations 
of relativistic stars and by Andersson and Th ylw e to de- 
scribe scattering from the Schwarzschild BHJ24| as well 
as to interpret the Schwarzschild BH glory [25]]. In this 
context, Andersson established, for the Schwarzschild BH 
of mass M, the existence of a family of "surface waves" 
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(each one associated with a Regge pole of the S-matrix) 
orbiting close to the unstable photon orbit at r = 3M (see 
also Ref. [l5[ for a more rigorous approach). Recently, 
from these "surface waves", we have been able to theo- 
retically and numerically construct the spectrum of the 
weakly damped complex frequencies of the Schwarzschild 
BH QNMs [TBI, [I?} and to interpret them as Breit-Wigner 
resonances. This has been achieved by obtaining ana- 
lytically the nonlinear dispersion relation as well as the 
damping of the "surface waves" propagating close to the 
photon sphere. Let us also note two related papers con- 
cerning analytical or numerical determinations of the 
Regge poles of the Schwarzschild BH [12, Hfil and that, 
this last year, Regge poles have also been used to under- 
stand the resonant aspects of the BTZ BH [l6| as well as 
to analyze some aspects of self- force calculations [27| ■ 

In the present paper, we extend the analysis devel- 
oped for the Schwarzschild BH to more general BHs and 
we establish, from Regge pole considerations, a precise 
connection between the existence of a photon sphere and 
the properties of the "surface waves" propagating close 
to it. More precisely, we consider a wide class of static 
spherically symmetric BHs of arbitrary dimension with 
a photon sphere, i.e., a hypersurface on which a mass- 
less particle can orbit the BH on unstable circular null 
geodesies. For more rigorous definitions of the photon 
sphere concept in static spherically symmetric spacetime, 
we refer to the article by Claudel, Virbhadra and Ellis 
28]. This class of BHs includes various spacetimes of 
physical interest such as Schwarzschild, Schwarzschild- 
Tangherlini and Reissner-Nordstrom BHs, the canonical 
acoustic BH or the Schwarzschild-de Sitter BH. For this 
class of BHs, we provide general analytical expressions 
beyond the leading-order terms for the Regge poles of the 
S'-matrix associated with a massless scalar field theory. 
These results permit us to obtain analytically the non- 
linear dispersion relation and the damping of the "sur- 
face waves" lying close to the photon sphere as well as, 
from Bohr-Sommcrfeld-type resonance conditions, the 
complex frequencies corresponding to the weakly damped 
QNMs. 

Our paper is organized as follows. In Sec. II, we display 
our general working assumptions and we justify them 
physically. We then explain how to construct the S'- 
matrix permitting us to analyze the resonant aspects of a 
scalar field theory defined on an asymptotically flat static 
spherically symmetric BH of arbitrary dimension with a 
photon sphere and we finally define its Regge poles as 
well as its complex quasinormal frequencies. In Sec. Ill, 
we provide a general analytical expression for the Regge 
poles. This is achieved by using and extending the WKB 
approach developed in the context of the determination 
of the QNMs by Schutz and Will [H and by Will and 



Iyer [30j,|3l| (see also Ref. [32J for general aspects of WKB 



theory and for particular aspects connected with eigen- 
value problems). Our result permits us to describe the 
Regge trajectories of a general asymptotically flat static 
spherically symmetric BH of arbitrary dimension with a 



photon sphere and to obtain, from semiclassical formulas, 
analytical expressions for the QNM complex frequencies. 
Our WKB analysis permits us moreover to show that (i) 
the dispersion relation of the nth "surface wave" is non- 
linear and depends on the index n and that (ii) the damp- 
ing of the nth "surface wave" is frequency dependent. 
In Sec. IV, we apply the general theory developed in 
Sec. Ill to particular BHs (Schwarzschild, Schwarzschild- 
Tangherlini, Reissner-Nordstrom and canonical acoustic 
BHs). In a brief conclusion, we consider some conse- 
quences of our work as well as possible extensions. In 
Appendix A, we establish the semiclassical connection 
between the Regge poles of a static spherically symmet- 
ric BH of arbitrary dimension with a photon sphere and 
the complex frequencies of its weakly damped QNMs. 
In Appendix B, we consider the particular case of the 
Schwarzschild-de Sitter BH. Indeed, even if such a grav- 
itational background is not asymptotically flat, the for- 
malism developed in Sees. II and III naturally applies to 
it. 

In this paper, we shall use units with h = c = G = 1. 



II. QUASINORMAL FREQUENCIES AND 
REGGE POLES OF STATIC SPHERICALLY 
SYMMETRIC BLACK HOLES: GENERAL 
THEORY 

We consider a static spherically symmetric spacetime 
of arbitrary dimension d > 4 with metric 



ds 2 



-f(r)dt 2 



dr 2 
W) 



r 2 da}_ 2 . 



(1) 



Here da 2 ,_ 2 denotes the line element on the unit sphere 
S d ~ 2 . On S d ~ 2 , we introduce the usual angular coordi- 
nates 6i € [0, 7r] with i = 1, . . . , d- 3 and ip € [0, 2tt}. We 
have 



do 2 d-2 




del. 



Yl sin 2 6i dip 2 

\i=l ) 



(2) 

Of course, a metric such as (JTJ) does not describe the most 
general static spherically symmetric spacetime but it will 
permit us to consider a wide class of BHs of physical 
interest. 

In Eq. ((TJ), we shall furthermore assume that f(r) is a 
function of the usual radial coordinate r with the follow- 
ing properties: 

• (i) There exists an interval I =]r^,+oo[c R with 
rh > such as f(r) > for r € J. 

• (ii) rh is a simple root of f(r), i.e., 

f(r h )=0 and f'(r h )?0, (3) 
and f(r) moreover satisfies 



lim f(r) = 1. 

r— >+oo 



(4) 
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(iii) There exists a value r c € I for which 



and 



/' (rc) - -/(r c ) = 

fc 



/"(r c ) - -j /(r c ) < 0. 



(5) 



(6) 



We shall now briefly discuss assumptions (i)-(iii) previ- 
ously introduced. Assumptions (i) and (ii) indicate that 
the spacetime considered is an asymptotically flat BH 
with an event horizon at r^, its exterior corresponding 
to r € J. Assumption (iii) implies the existence of a 
photon sphere which is the support of unstable circular 
null geodesies (see below for more details). It should be 
noted that, as a consequence of (i) and (ii), the tortoise 
coordinate r* = r*(r) defined for r € I by the relation 
dr t /dr = l/f(r) and the condition r + (r c ) = provides a 
bijection = r*(r) from I to ] — oo, +oo[. 

Let us consider a free-falling massless particle orbiting 
the BH. Without loss of generality, we can consider that 
its motion lies on the equatorial hyperplanc defined by 
6i = 7r/2 for i = 1, . . . , d — 3. Because it moves along a 
null geodesic, we have [cf. Eqs. ([I) and (|2j] 



dt 
da 



f(r) 



dr 
da 



dip 
da 



(7) 



where a is an affine parameter and, of course, there exist 
two integrals of motion respectively associated with the 
Killing vectors d/dt and d/dip and given by 



f(r) 



dt 
da 



E, 



dy 

da 



= L. 



(8a) 
(8b) 



Here E and L denote respectively the energy and the 
angular momentum of the massless particle. Inserting 
Eqs. (JHaJ) and j8b| into Q, we obtain 



where 



Ta) 



V cS (r) = ^/(r) 



E z 



(9) 



(10) 



From these last two equations and from assumption (iii) , 
one can easily remark that the massless particle can orbit 
the BH on an unstable circular geodesic defined by r = 
r c . Indeed, we have in particular 



dr 



V eS (r) 







and 



£ (/"(rc) " £/(r„) 



(11a) 



<0. (lib) 



On this orbit, the massless particle takes the time 

2nr c 



T 



(12) 



to circle the BH. This result can be obtained by integrat- 
ing Eq. ©. 

The wave equation for a massless scalar field propagat- 
ing on a general gravitational background is given by 



1 



-A, 



-9 



-gg^d^)=Q. (13) 



If the spacetime metric is given by ([T]), after separation of 
variables and the introduction of the radial partial wave 
functions 3>g(r) with £ = 0,1,2,..., this wave equation 
reduces to the Regge- Wheeler equation 



d 2 t> e 
dri 



+ [uj 2 - Vt(r)] §i = 0. 



(14) 



[Here we have assumed a harmonic time dependence 
exp(— iujt) for the massless scalar field.] In Eq. (flT)) . Vt{r) 
is the Regge- Wheeler potential given by 



Ve(r) = f(r) 



£(£ + d~3) (d-2)(d-4) 



4r 2 
d-2 
2r 



f(r) 



fir) 



(15) 



It should be noted that 



• lim r _> rfl Vi (r) = and lim r _ ) . +00 Vi (r) = and 
therefore the solutions of the radial equation (TT4")) 
have a behavior in exp(±iu;r*) at the horizon and 
at infinity. 

• For £ 3> 1, Ve(r) has a local maximum at r = r c 
because, in this limit, Vg{r) and V c $(r) are similar. 

• For any finite value of £, the local maximum of Vi (r) 
is close to r = r c . 

For a given angular momentum index £, the 5-matrix 
element Si(uj) is defined by seeking the solution of the 
Regge- Wheeler equation (TI"4"|) which has a purely ingoing 
behavior at the event horizon r = r^, i.e., which satisfies 



$t(r) - T t (u)e- 

r. M — ¥ — oo 

and which, at spatial infinity r — 
asymptotic behavior of the form 



(16) 

-oo, presents an 



$ £ (r) - < 

r* — f+oo 



-iwr,+i(^+^2)f-if 

S l {uj)e +lbjr '< l+ ^ L )^ +l T . (17) 



We recall that the 5*-matrix permits us to analyze the res- 
onant aspects of the considered BH as well as to construct 
the form factor describing the scattering of a monochro- 
matic scalar wave (see Appendix A). 
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To describe semiclassically resonance phenomena, the 
dual structure of the 5*-matrix plays a crucial role. In- 
deed, the S'-matrix is a function of both the frequency w 
and the angular momentum index I. It can be analyti- 
cally extended into the complex w-plane as well as into 
the complex A-plane (CAM plane) with A = £+(d— 3)/2. 
Prom now on, we shall denote by 5 , A-(d-3)/2( a; ) this dou- 
ble analytical extension. For I = A — (d — 3)/2 G N, 
the simple poles lying in the fourth quadrant of the 
complex w-plane [let us recall that they are also simple 
poles of Tx-(d-3)/2( UJ )} are the complex frequencies of the 
QNMs. These modes are therefore solutions of the radial 
wave equation (|14[) which are purely outgoing at infin- 
ity and purely ingoing at the horizon. We shall denote 
by Loe n = uj^n — iT£ n /2 where n € N* the quasinormal 

frequencies. We recall that uj^ > and T[ n > repre- 
sent respectively the frequency of the oscillation and the 
damping corresponding to the associated QNM. We as- 
sume that, in the immediate neighborhood of LU£ n , Si(uj) 
has the Breit-Wigner form, i.e., 



Se(w) oc 



T^n/2 



,(°) 



(18) 



For a given value u > of the frequency, the simple poles 
lying in the first quadrant of the complex A-plane are the 
so-called Regge poles. It should be noted that they are 
also poles of T X -(d-3)/2{ UJ ) an d therefore the associated 
modes (Regge modes) are purely outgoing at infinity and 
purely ingoing at the horizon. We shall denote the Regge 
poles by A„(w), the index n = 1, 2, . . . permitting us to 
distinguish each pole. 

The structure of the S'-matrix in the complex A-plane 
allows us, by using integration contour deformations, 
Cauchy's theorem and asymptotic analysis, to provide 
a semiclassical description of scattering (see Appendix A 
for more precisions). The curves traced out in the CAM 
plane by the Regge poles as a function of the frequency 
bj are the so-called Regge trajectories. They permit us 
to interpret Regge poles in terms of "surface waves" (see 
Appendix A): ReA n (w) provides the dispersion relation 
for the nth "surface wave" while ImA n (w) corresponds 
to its damping. Furthermore, from the Regge trajecto- 
ries, we can semiclassically construct the resonance spec- 
trum [see, in Appendix A, formulas (|A9|) . (jAllI) and 
(IA12[) ]. The semiclassical formula (a Bohr-Sommerfeld- 
type quantization condition) 



Re A, 



d- 3 



G N 



(19) 



provides the location of the excitation frequencies ui^ of 
the resonances generated by nth "surface wave" , while 
a second semiclassical formula gives the widths of these 
resonances 



Im X n (oj) [d/ du Re X n (oj)] 



[d/dujReX n (uj)} 2 + [d/duj Im A„(w)] 2 



(20) 



It should be moreover noted that this formula re- 
duces, in the frequency range where the condition 
\d/dujRe X n (ui)\ » \d/duj Im A ra (w)| is satisfied, to 



Tin _ Im X n (u>) 
2 d/duiRe A„(cj) 



jo) 



(21) 



III. WKB APPROXIMATIONS FOR THE 
REGGE POLES AND SEMICLASSICAL 
EXPRESSIONS OF THE COMPLEX 
QUASINORMAL FREQUENCIES 



In general, it is not possible to solve exactly the Regge- 
Wheeler equation (IT4l) and therefore we can obtain only 
analytical approximations for the Regge poles and for 
the complex quasinormal frequencies. For example, the 
WKB approach developed in the general context of eigen- 
value problems (for more details see Ref. [HI]) has been 
adapted for the determination of the Schwarzschild BH 
QNMs by Schutz and Will and by Will and Iyer 
[30l HH and, in Ref. [l7j], for the determination of the 
Schwarzschild BH Regge poles. It can be extended to 
the more general case considered in this paper. By using 
third-order WKB approximations [30l . l31j for the Regge 
modes of Eq. (IT41 , we find that the Regge poles A are the 
complex solutions of the equation 



T/ (A) + [-2F (2) (A)] 1/2 A(A,n) 



-ia(n) [-2t/ (2) (A)] 1/2 [1 + n(A, n)] (22) 



with u > and n — 1, 2, 3, Here 



A(A,n) 



-2^ (2) (A)] 
288 {V^(X)J 



1/2 



8y (2) (A)l4 [} 



(7 + 60a(n) 2 



(23a) 



5 



and 



[-2F (2) (A)] 
5 /V (3) (A)\ 



6912 V^ (2) (A), 



(77+ 188 a(r 



1 / [V^(\)YV^(X) 



K, (4) (A) 



+ 



384 
1 

'2304 ^(A), 

1 ^o (3) (A)F (5) (A) 
288 



(51 + 100a(n) 2 ) 



(67 + 68a(n) 



'l9 + 28a(n) 



1 /V (6) (A)* 



(5 + 4a(n) 2 ) 



(23b) 



288 V^ (2) (A) / 
In Eqs. (f2"2"j) and (f2U)l . we have introduced the notations 
a(n) = n- 1/2 (24) 

and, for p £ N, 



K, (p) (A) = -^ A _ (d _ 3)/2 (r*) 



dr 



r« = (r,) 



(25) 



with (r*) which denotes the maximum of the function 
VA_(d_ 3 )/ 2 (r«). 



In order to solve Eq. (|22j) . we need to express the 
asymptotic expansions for |A| — > +oo of Vq(X) and 



2y (2) (A) 



1/2 



and of the various ratios appearing in 

Eqs. (|23al) and J23b|. In order to simplify the results we 
have introduced the notations 



/W = /W( rc ) 



(p)i 



and 



,,, = -\ L/ c 2rH/. 



(2) 



(26) 



(27) 



It is worth noting that the rj c parameter is directly linked 
to the second derivative (jllb[) of the effective potential 
(|T0|) taken at r = r c . As a consequence, it represents a 
kind of measure of the instability of the circular orbits 
lying on the photon sphere. In fact, it can be expressed 
in terms of the Lyapunov exponent A c corresponding to 
these orbits introduced in Ref. [llj and which is the in- 
verse of the instability time scale associated with them: 
we have 



Vc 



:|A C 



(28) 



We will say no more about this connection because, as al- 
ready mentioned in Sec. I, we intend to go beyond purely 
geometrical considerations in our analysis of the resonant 
behavior of BHs. 

After a tedious calculation, we obtain 



and 



as well as 



-2V (2) (A) 



2r) c f c 



2/ 2 [(d-3) 2 + (d-2)(d-4)/ c 



+r c 2 / c /j 2 )[(d - 2)(d + 8)/ c - 2(d - 3) 2 ] + d(d - 2)rlf 2 J^ 



+rt (/ c (2) ) 2 [(l/2)(d-3) 2 -(d 2 -4)/ c 



a|^+oo yA- 3 



y (2) (A) 



fc 



2rf c r 2 c 

+4 U(fpy+fjw 



16/ 2 -16r c 2 / c /i 2 )+4r c a / c /, 



f(3) 



° T2 

A|^+oo V X z 



(29a) 



(29b) 



(30a) 



A|^+oo \X Z 



(30b) 



6 



^o (3) (a)]V> ( a) ____ rm(^y 



8r?| 

+ri (4 (/J 9 ))' + /«/«) 



16/ c 2 -16r c 2 / c / c (2) +4r3/ c /( 3 ) 



° U2 I • 

a|^+oo VA- 2 / 



(30c) 



ltf 3) (A)V < B >(A) _ r c 2 / 3 / c (3) 



y (2) (A) 



4?7 4 



-Wf c fP + WrJ c fW 



(l5/< 2 >/< 3 > + /c/ c (5) 



° \ v 

|A|->+oo \\ z 



(30d) 



K) (6) (A) / 2 

Vq (2) (A) 2f 7 2 '< 



-272/ c 3 + 408r 2 / 2 / c (2) - 88r 3 / 2 / c (3) 



+rtfc (38/ c /i 4 > - 204 (frf ) + r c 5 / c (l04/( 2 )/j 3 ) + 18/ c/ ( 5 >) 



+r° 34 



(/i 2 )) 3 + 15/ c (/J')) 2 + 26 /c /i 2 )/i 4 ) + / 2 /i 6 ) 



(30e) 



I 

We can now solve Eq. (|2"2"|) by assuming |A| ^ 1 as well A„(w) with n = 1, 2, 3, . . . given by the approximation 
as Re A 3> ImA. We obtain for the solutions a family 



where 



and 



A n (cj) 



■ u} 2 + a n + 2rj 2 a(n) 2 e n (oj) 



1/2 



+ irj c a(n) [1 + e„(w)] 



1 = ~ nkf c { 288/ ' [(d2 " 2d ~ 1)/c " (d ~ 3)2] 

+ 144r 2 / c /i 2 ) [2{d - 3) 2 - (2d 2 -Ad- 3)/ c ] 



-72rf/ c 2 /i a) - 18r 4 
+36r c 5 / c /i 2) / c (3) + ^ 



4(d - 3) 2 (/| 2 )) 2 - 4(d - 3)(d + l)/ c (/ c (2) ) 2 + / 2 / c (4) 
36 (/ c (2) ) 3 -7/ c (/i 3 )) 2 + 9 /c /i 2 )/i 4 ) 



+ «W 2 |f { 24 / C /i 3) + 6-c/ c /i 4 » - 12r 2 /( 2 >/i 3 > + r 3 ( 5 (/i 3 ^ - 3/j 2 >/i 4 >) 



(31) 



(32) 



{r 2 /f c )u 2 + a n + r]%a(n) 2 



(33) 



with 



b n = 



442368ry, 

-9216(3d 2 - 6d - 2)^/ c 5 / c (3) + 6912r 4 /, 

-192r c 6 / 2 
-5/ c 3 / c (6) " 



3456r c 5 / c 4 



+ 12r c 8 /c 



+48r c 9 / c 2 



+ 12r, 



^ {-55296(d- 3) 2 / c 5 (l + f c ) + 27648r c 2 / c 4 / c (2) [5(d- 3) 2 + 2(d 2 - 12d + 23)/ c ] 

20(d- 3) 2 (/ c (2) ) 2 +2(2d 2 + 36d- 95)/ c (/ c (2) ) 2 
6(2d 2 -4d+l)/i 2 )/( 3 )+5/ c /i 5 )" 
360(d - 3) 2 (/ c (2) ) 3 + 144(2d 2 + 6d- 25)/ c (/J 2 ))' + 423/ 2 /( 2 '/ c (4) - 200/ 2 (/J 3 ) 
96r c 7 / 3 36(6d 2 - 12d + 17) (/ c (2) ) ' / c (3) - 257 fj^f^ + 270 / c / c (2) / c (5) 
-1440(d- 3) 2 (/P) 4 + 288(7d 2 + 6d - 55)/ c (/ c (2) ) 4 + 2376/ 2 (/j 2 >) 2 /<*> 
+67/ 3 (/i 4 )) 2 + 152/ c 3 /i 3 )/f) - 2744/ 2 /i 2 ) (/j 3 >) 2 - 120/ c 3 /?>/j 6) 
12(6d 2 - 12d + 59) (/i 2 )) 3 /i 3) + 153/c (/ c (3) ) 3 - 514/ c /i 2 >/( 3 )/j 4 > + 270/ c [f^ /j 5 > 
144(d- 3) 2 (/i 2 )) 5 - 288(d 2 - 7)/ c (/j 2 >) 5 - 67/ c 3 /j 2 ) (frf - 152 /c 3 /c (2) / c (3) / c (5) 
+344/ 2 (/i 2 >) 2 (/i 3 )) 2 + 60/ 3 (/P) 2 /(«) + 108/ 2 (/P) 3 /W + 153/ 3 (/j 3 >) 2 /« 



-24r 41 / c / c (2) 
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252 (/( 2 >)' /f) + 153/ c (/i 3 )) - 257/ c /i 2 )/( 3 )/i 4 ) + 90/ c (/( 2 >) 
864 (/i 2 )) 6 + 385/ 2 (/F) 4 - 1512/ c (/| 2 >) 4 /M - 91 8/ 2 /( 2 > (fP)" fW 
+201/ 2 (/P) 2 (/W) 2 + 456/ 2 (/i 2 )) 2 /i 3 >/( 5 > + 1368/ c (/( 2 >) 3 (/W) 3 - 120/ 2 (f^f 
• 11059^10 { 9216 /c/ c (3) + 13824r c / 4 /j 4 > + 3456r c 2 / 3 [-2 /c ( 2 )/( 3 > + / c / c (5) 



+192r 3 / 3 
+ 12r c 5 /, 



5 f2 
cJc 



72 (/f)) _99/(2)/W+/ c /f) 
648 



288r 4 / 2 



12 (/i 2) ) 2 /< 3) - 29/ c /i 3 >/W + 18/ c /i 2 )/i 5 ) 



144r c 6 / c 



-12r c 7 / c 



(/ c (2) ) /i 4) + 17/c (/ c (4) ) +56/ c /i 3 )/i 5 )- 1032/i 2 )(/i 3 )) -24/ c /( 2 )/i 6 ) 
28 (/P) 3 /J 3 > + 25/ c (/f)) 3 - 58/ c /( 2 )/i 3 ) /c ( 4 ) + 18/ c (/j 2 >) 2 /(« 

/i 2) ) 3 /i 4) + 75/ c (/i 3 )) 2 /i 4 ) - 17/ c /( 2 ) (/i 4 )) 2 - 56f c fPfPf& 

/i 2) ) 2 (/c (3) ) 2 + 12 /c(/ c ^) 2 /^ 



-36 
+ 168 



-72r c 8 / c (2) 



12 



(/i 2 )) 3 /i 3 ) + 25 /c (/i 3 )) 3 - 29 /c /i 2 )/i 3 )/i 4 ) + 6/c (/ c (2) ) 2 / c (5) 



235/ c (/i 3 )) 4 - 216 (/i 2 )) 4 /« - 450/ c /i 2 ) (/j 3 >) 2 /j 4 > + 51/ c (/j 2 >) 2 (/W 
+168/ c (/P) 2 /i 3 )/i 5 ) + 360 (/i 2 >) 3 (Ji 3 ') 2 - 24/ c (/P) 3 /J 6 ) 



(34) 



Equation (131 [) is the main result of our paper. As we shall see below (see also Appendix A), it also provides 
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expressions for the dispersion relation and the damping 
of the "surface waves" lying on (close to) the photon 
sphere of the considered BH. 



It is moreover possible to simplify (|3Tj) and to obtain 
a high frequency approximation for the Regge poles. We 
have 



J 



X n {uj) 



(2r c ///Q lo 



iri c a{n) 
I — 



O 

oj— >-+oo 



(35) 



It is important to understand the difference between the 
approximation pip and its much more elegant version 
(|35p . The approximation (I35[) is meaningful as a 1/u ex- 
pansion with lo — > +00. By contrast, the approximation 
(I3ip remains valid in a large range of frequencies thanks 
to WKB theory. Indeed, in order to establish Eq. ([22]) 
we have considered as a perturbation parameter of the 
WKB method the distance between the turning points of 
Eq. CEI]) [i.e., the roots of Qe(r*) = lo 2 - V^(r*]land the 
location of the peak of Qt{r*) (see also Refs. 29 32f ) in- 
stead of 1/uj. Of course, in order to solve Eq. (|2"21 and to 
obtain the expression pip , we have furthermore assumed 
that |A| ^> 1 as well as Re A ^> ImA and, as a conse- 
quence, we cannot expect from pip a very high accuracy 
for very low frequencies (see also Refs. [l2[ and [l?} for 
related numerical studies in the case of the Schwarzschild 
BH). 

It is furthermore interesting to provide an interpreta- 
tion of the previous results in terms of "surface waves" 
as it is customary in the CAM approach. As we have 
noted in Appendix A, the resonant part of the form fac- 
tor is a superposition of terms like e~xjp[i(\ n (oj)6i— cot)] (in 
this paragraph we take into account the harmonic time 
dependence exp[— icot]). By inserting the leading-order 
terms of Eq. p5p into these wavelike contributions, we 
can easily note that the contribution of the nth "surface 
wave" reduces to 



exp[— r] c a(n)0i] exp 



(rch 



f c v6 



■ Lot 



(36) 



The second term describes the propagation of this "sur- 
face wave" near the photon sphere at r — r c . Indeed, 
it circles the BH in time T" = 27r(r c / V7c) which is ex- 
actly the time (|T2"j) needed for a massless particle to orbit 
the BH on an unstable circular null geodesic. The first 
term corresponds to an exponential decay of this "surface 
wave" due to continual reradiation of energy. It is very 
interesting to rewrite Eq. p6p in the form 

(37) 



cxp[— k[[(uj)C] exp[i(k' n (uj)£ — uxb)] 



with C = r c 6\ which denotes the arc length taken on 
the photon sphere. Now, k' n {uj) = Lo/\ffc represents the 
wavenumber of the nth "surface wave" or, in other terms, 
its dispersion relation, while fc"(cj) = n c a(n)/r c is its 
damping constant. 

To the leading order, the dispersion relation is linear 
and independent of the index of the "surface wave" while 



the damping constant depends only on the index n. Of 
course, if we go beyond the leading-order terms, the dis- 
persion relation reads [see Eq. p~Tl) ] 



LO 



-,1/2 



a(n) 2 e„(w) 



(38a) 



which implies [see Eq. p5l) ] 



*4M = 



o ^ 



(38b) 

and is clearly nonlinear as well as dependent on the in- 
dex n. This result could have important consequences 
in strong gravitational lensing (see also Ref. [ljj and the 
discussion in Appendix B.l of the present paper). With 
such a potential application in mind, it is worth noting 
that, in Eqs. P"T|) and p"5|) . lo denotes the frequency of the 
scalar photon observed by a static observer at infinity. If 
we consider the gravitational redshift of this photon and 
introduce its frequency Co — lo/ \J~fc measured by a static 
observer lying on the photon sphere, the dispersion rela- 
tion (l38l) reads 



and we have 



Co 2 



^- + 2 



2r? Co 



° I" 

LO 6 



1/2 



(39a) 



(39b) 



When a n > 0, Eq. p^| provides a superluminal dis- 
persion relation. Indeed, it leads to a group velocity 
v g (Co) = dCb/dk' n {Co) > 1. The condition a n > is 
satisfied Vn e N for spins 1 and 2 propagating on the 
Schwarzschild BH [see Eq. (13) of Ref. H3]. It is also 
satisfied for the massless scalar field when 11 / 0. As a 
consequence, it cannot be considered as an exotic phys- 
ical condition. It even seems to be true in most cases 
encountered (see the examples considered in Sec. IV). At 
first sight, this result may seem a little bit puzzling. But 
it is important to note that superluminal dispersion rela- 
tions for the "surface waves" lying on the photon sphere 
do not necessarily lead to a violation of the relativistic 
principle of causality. Indeed, for the transfer of infor- 
mation between a source and a receptor located outside 
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a BH, various "channels" are involved. Of course, there 
are channels associated with diffraction by the BH pho- 
ton sphere [for a source and a receptor at infinity, they 
correspond to the sum over the Regge poles in Eq. (|A8|) ] 
but there are also channels associated with geometrical 
rays [for a source and a receptor at infinity, they come 
from the background integral (IA2[) over the contour V of 
Fig. [T] after asymptotic evaluation]. In order to study 
causality, it would be necessary to carefully take into ac- 
count interferences between all the monochromatic spec- 
tral components of the signal carrying the information for 
all the channels involved. We believe that such a study 
would not show a violation of causality. The situation 



encountered here is similar to that discussed by many 
authors working in electromagnetism of dispersive media 
(see, e.g., Ref. [331 and references therein). In such a con- 
text, it has been observed that a group velocity greater 
than the light velocity does not violate causality because 
it is not the velocity of information transmission or the 
energy velocity. 

Finally, the WKB result (|3T1) and the associated Regge 
trajectories permit us to derive, from the semiclassical 
formulas (fT9"f and (PHI) , useful analytical expressions for 
the QNM complex frequencies. Indeed, by inserting ([31]) 
into Eqs. (fT9|) and (|2Tj) . we obtain the large £ behaviors 



(o) _ v%. 



d-3\ a n a n 



-Soo \£ 3 



- in y/Yc / \ 

— = Vc a{n) 

Z Tr 



with £ € N and n = 1,2,... Here 



2£ 2 i-^+oo V £ 3 



c n —- a n -f- 26 n . 



(40a) 
(40b) 

(41) 



r 



It should be noted that the leading-order terms of 
Eq. (j40l) have already been obtained in Refs. ficl - ll^ . 
The higher-order terms in l/£ and l/£ 2 are new (see, 
however, Sec. 5.2 of Ref. [Hj]). They are directly linked 
to the nonlinear behavior of both the dispersion relation 
and the damping of the "surface waves" lying on (close 
to) the photon sphere. 



IV. APPLICATIONS 

In this section we shall apply the previous formalism 
to various spacetimes of physical interest. 



A. Schwarzschild-Tangherlini black holes 

For Schwarzschild-Tangherlini BHs [33], the function 
f(r) reads 



The unstable circular null geodesies are located at 



f(r) = 1 - 



d-3 



with 



-d-3 



16ttM 



(d - 2)A d - 2 



and Ad- 



2 ^-i)/2 

r[(d-i)/2] 



(42) 



(43) 



Here M is the mass of the BH, Ad-2 the area of the unit 
sphere S d ~ 2 and the event horizon is located at = ro- 



r c = r 



d-1 



l/(d-3) 



(44a) 



(44b) 



and the associated rj c parameter is given by 
f] c = yd — 3. 

1. d=4- The Schwarzs child black hole 



For d = 4, the Schwarzschild-Tangherlini solution is 
nothing but the ordinary Schwarzschild BH and we have 



as well as 



r c = -r = 3M, 
Vc = 1, 



29 5 . , 2 

1 ainYi 

216 18 w 

1357 305 , , 2 
-a(n) , 



15552 3888 
313 235 



7776 1944 



a(n) 2 , 



27M 2 w 2 + a„ + a{nf 



(45a) 
(45b) 



(46a) 
(46b) 
(46c) 
(46d) 
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The WKB approximation (|3Tj) for the Regge poles leads 
to 

A„(w) w [27M 2 w 2 + a n + 2a(n) 2 e„(w)] 1/2 

+ia(n)[l + e„(w)] (47) 

and their high frequency behavior (1351) provides 

V3a r , 



A„(w) 



3a/3Mc 



+ia(n) 



1 



18Mw 



27Af 2 w 2 



P, 73 M 48 ) 



Formulas (14T1) and (I48p are in agreement with the results 
obtained in Ref. [17j . 

The resonance excitation frequencies and the damping 
of the QNMs given by the general formulas (|4T)|) reduce 
to 



1 



3V3M 

a(n) 
3V3M 



1 



2£ 2 



2t At 2 t^+oo \ P 



1 



O 



(49a) 
(49b) 



Formulas (49j are in agreement with the results obtained 
in Ref. [31j(see also Ref. [13]). 



2. The five-dimensional Schwarzschild-Tangherlini black 
hole 



For <i = 5, we have 



as well as 



r c = V2r , 
Ve = >/2, 

Q 

a(n) 2 , 



a(nf 



16 4 
115 31 

512 ~ 128 

5 17 , si 

I an , 

48 64 y ' 1 



£n(w) 



2r 2 w 2 + a„ + 2a(n) 2 



(50a) 
(50b) 

(51a) 
(51b) 
(51c) 
(51d) 



The WKB approximation (j3"Tj) for the Regge poles leads 
to 

X n (oj) [2r 2 w 2 + a„ + 4a(n) 2 e„(u;)] 1/2 

+io(n)[l + e»(a;)] (52) 

and their high frequency behavior (|35|) provides 



A„(cj) = 



72 



2\/2r c . 



-i\2a{n) 



n 




+ O 






2r 2 w 2 _ 







The resonance excitation frequencies and the damping 
of the QNMs given by the general formulas (|4l7| reduce 
to 



1 



a(n) 



1 + 2P + £->+oo ( £3 



(54a) 
(54b) 



3. The six-dimensional Schwarzschild-Tangherlini black 
hole 



For d = 6, we have 



1 2 



1/3 



ro, 



as well as 



a " = ~60 + 5 a(n) ' 
1411 91 2 

0„ = -r-nr-r - — — a{n) , 



3600 180 
481 7 



1800 18 



a(n) 2 , 



£n(w) 



(55a) 
(55b) 

(56a) 
(56b) 
(56c) 
(56d) 



(5/3)r 2 u; 2 + a„ + 3a(n) 2 ' 
The WKB approximation (|3"Tj) for the Regge poles leads 



to 



A„(w) 



r c ui + a n + 6a(n) e n {tu) 



1/2 



+iV3a(n) [1 + e„(w)] (57) 
and their high frequency behavior (|35p provides 



A„(w) = 
+iv / 3 a(n) 



1/2 



(5/3)7 



2(5/3) 1 / 2 r c w 



O 



oj— y+oo \ u 



(58) 



The resonance excitation frequencies and the damping 
of the QNMs given by the general formulas (|40"1) reduce 
to 



,(°) 



~2 



3\ 1/2 M 



2^ 4£ 2 e^+oo W 3 



3\ 1/2 fl 



a(n) 



5 I \ r, 

T 2€ 2 £->+oo V P 



(59a) 



(59b) 
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4- Leading- order terms for the Schwarzschild-Tangherlim 
black hole of arbitrary dimension 

It should be noted that, in the d-dimensional case, we 
can easily derive the leading-order terms of the Regge 
poles and quasinormal complex frequencies. We have 



I 



, o r c uj + Wd-3a{n) + O — 
a — 3 uj->+oo V to 



(60) 



and 



,(°) - 



d- 3 m ( d _-3 



d - 1 \r, 



O 

-S- + 00 



2 Vd~^T \r c J y ' \l 



(61a) 
(61b) 



Formulas (I61|) are in agreement with the results obtained 
in Refs. [35j and ficj . 



B. The Reissner-Nordstrom black hole 

For the d-dimensional Reissner-Nordstrom BH 34] , the 
function f(r) reads 



/W = 1 -(7)"" + ^ 



(62) 



with (see Ref. [36| or Appendix A of Ref. [37] ) 

87^ 



16ttM 

(d - 2)Ai_ 2 



2/i and 9~ 



(d-2)(d-3)' 
(63) 

Here M is the mass of the BH and (7 is its charge. We 
furthermore assume that /i 2 > 8 2 . It should be noted 
that here we consider electromagnetism in the Heaviside 
system of units. For this background there are two hori- 
zons, the so-called inner and outer horizons. They are 



respectively located at 

r + = (/1 + VA^— ¥) 



l/(d-3) 
l/(d-3) 



(64a) 
(64b) 



We are only interested in the outer horizon with radius 
at Th — r + because we have /(r) > for r G]r/,, +oo[ [see 
assumption (i) of Sec. II]. 

The unstable circular null geodesies are located at r c > 
r h given by 

l/(d-3) 



(d - 1)M + v /[(^-l)M] 2 -4(d-2)0 2 



and we have for the associated parameter r\ c 



(65a) 



V ?*c 

For = we recover the Schwarzschild-Tangherlini BH 
results. It should be noted that it is also possible to 
express r\ c in two other equivalent forms since, as a con- 
sequence of Eq. ([5]), the parameters r c , fx and are linked 
by 

(rt 3 Y - (d- 1)K?~ 3 + (d - 2)6» 2 = 0. (66) 



1. The four- dimensional Reissner-Nordstrom black hole 



For d = 4, we have 

r c = - (3M + ^9Af 2 - 8Q 2 ) , (67a) 



1 2Q2 



(67b) 



as well as 



K = 



(i-Q 2 A 2 ) 2 

216(1 -2Q 2 /r 2 ) 2 

(i-Q 2 A c 2 ) 2 

15552(1 - 2Q 2 /r 2 ) 5 
1 - Q 2 /r c 2 



1 - 2Q 2 /r 2 ) 5 

(i-Q 2 A 2 ) 2 

" 7776(1 - 2Q 2 /r2)5 

1 - Q 2 /r c 2 
1944(1 - 2Q 2 /r 2 ) 5 



-29- 

1357- 
305 



21 



36(l-2Q 2 /r 2 ) 2 



10- 



62Q 2 124Q 4 



a(n) 2 , 



12478Q 2 42333Q 4 64324Q 6 43588Q 



313 - 
235 - 



V c T c T c T c 

3943Q 2 21335Q 4 | 56357Q 6 69544Q 8 | 34772Q 

T c T c T c T c 

3118Q 2 11229Q 4 18820Q 6 18820Q 



,10 



(68a) 



a(n) 2 , (68b) 



iyi2 ^4: 'y'O 

c c c c 

2645Q 2 11929Q 4 28315Q 6 37592Q 8 18796Q 10 



,,10 



a(n) 2 , 



3c) 



12 



and 



bn 

3r c 2 /(l - Q 2 /r c 2 ) w 2 + a n + (1 - 2Q 2 /r c 2 )a(n) 2 ' 



(68d) 



Here we have noted Q 2 = Anq 2 . This permits us to compare our results with those for which electromagnetism is 
expressed in the Gaussian system of units. The WKB approximation ([3T]) for the Regge poles leads to 



A„(w) 



3r 2 



1 - Q 2 /r; 



:LJ 2 + On + 2 1 



2Q 1 



a(n) 2 e n (u>) 



1/2 



/ 2Q 2 
+ i\l 1 — a(n)[l + e n (u;)] 



(69) 



and their high frequency behavior (|35|) provides 



complex frequencies. We have 



X n {uj) 



V3i 



y/l - Q 2 /r 2 2V3r c /^l-Q 2 /r 2 u J 



■ 1 2 ^ ( \ 



o : 1 



i 



3r 2 /(l-Q 2 /r 2 )w 2 



(70) 



A„(w) = 



d- 1 



d- 3 



1 - 



-Hv d — 3 A / 1 ^r^— a(n) 



2(d-3) 



and 



o - 



(72) 



The resonance excitation frequencies and the damping 
of the QNMs given by the general formulas (|4"U|) reduce 
to 



Q 2 



r 2 \ V3r c 



e + 



21 U 2 <-v+oo I £ 3 



= 4 1 



2Q 2 /Q 2 ( 1 



Va/3) 



(71a) 



a(n) 



(71b) 



,(o) 



d-3 



'i- 



i 



2(d-3) V r 



o 

->+oc 

d- 3 



x/d - T 



d-3 



(73a) 



2 



2(d-3) 



(d-2)6 2 

2(d-3) 



(73b) 



It should be noted that for 9 = 0, from ([72} and ([75]). 
we recover the Schwarzschild-Tanghcrlini BH results of 
Sec. IV.A.4. 



The leading-order terms of (|71|) arc in agreement with 
formulas (37) and (38) of Ref. Q. Furthermore, it should 
be noted that, for Q = 0, all the results obtained for the 
Regge poles and the complex quasinormal frequencies of 
the four-dimensional Reissner-Nordstrom BH reduce to 
the Schwarzschild BH results of Sec. IV.A.l. 



2. Leading-order terms for the Reissner-Nordstrom black 
hole of arbitrary dimension 



In the d-dimensional case, we can easily derive the 
leading-order terms of the Regge poles and quasinormal 



C. The canonical acoustic black hole 

For the canonical acoustic BH (for more details see 
Sec. 8 of Ref. .38]), the function f(r) reads 



(74) 



The sonic event horizon is located at rn = ro, the radius 
of the unstable circular null geodesies is given by 

r c = 3 1/4 r (75a) 

and we have for the corresponding rj c parameter 

Vc = 2. (75b) 
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Furthermore, we have 

1 20 



an= io8 + y a(n)2 ' 

85 215 2 

On = an , 

1944 243 V ' ' 



19 110 

243 + 243 



a(n) 2 , 



(3/2)r 2 cj 2 + a n + 4a(n) 2 



(76a) 
(76b) 
(76c) 
(76d) 



The WKB approximation (|3"Tj) for the Regge poles leads 
to 



X n (uj) 



r c ui + a n + 8a(n) e n (tu) 



1/2 



+2ia(n)[l + e n {uj)} (77) 
and their high frequency behavior (|35|) provides 

X n (uj) 

+2ia(n) 



V2 V2r c oj 
1 



(3/2)7 



O 

uj— ¥-\-oo \ UJ 



(78) 



The resonance excitation frequencies and the damping 
of the QNMs given by the general formulas (14TJ1) reduce 
to 



and in the wave/field theory context, the appealing and 
intuitive interpretation of Schwarzschild BH QNMs sug- 
gested by Goebel in 1972 i.e., that they could be 
interpreted in terms of gravitational waves in spiral or- 
bits close to the unstable circular photon/graviton orbit 
at r = 3M which decay by radiating away energy. 

In the present paper, we have greatly extended the 
approach initiated in Ref. (l7| . More precisely, we have 
provided general analytical formulas for the Regge poles 
of the S'-matrix associated with a massless scalar field 
theory defined on a static spherically symmetric BH of 
arbitrary dimension with a photon sphere. This has been 
achieved by using third-order WKB approximations to 
solve the associated radial wave equation and by empha- 
sizing more particularly the role of the photon sphere or, 
in other words, of the unstable circular null geodesies on 
which a massless particle can orbit the BH. But, it is 
important to note that our results are only a first step 
to understand, in a semiclassical framework, various as- 
pects of static spherically symmetric BH physics such 
as wave scattering, gravitational lensing, Hawking radi- 
ation... Moreover, we think that it would now be very 
interesting to introduce the CAM method in a less sym- 
metric situation, e.g., to study Kerr and Kerr- Newman 
BHs (see Ref. [26| for a first numerical step in this direc- 
tion) and to analyze the splitting of their complex quasi- 
normal frequencies or to understand, from a semiclassical 
point of view, the superradiance phenomenon. We leave 
such a study for the future. 



,(°) 



y/2 



Tin 2\f2 



1 + ^-+ o (— 

2£ 2 1^+00 W 3 



(79a) 
(79b) 



The leading-order terms of (|79p are in agreement with 
formula (53) of Ref. In Ref. [13, Dolan and Ottewill 
have obtained for n = the expansions of u;!2 and I^„/2 
up to order l/(£+ 1/2) 4 . Our results (1791 are consistent 
with their Eq. (78). 



V. CONCLUSION 

As noted by Chandrasekhar in the mid-1970s, BH per- 
turbation theory can be formulated as a resonant scat- 
tering problem. As a consequence, all the techniques 
developed in the framework of scattering theory can be 
naturally introduced in the context of BH physics. One 
of the central concepts of scattering theory is the concept 
of a Regge pole which permits one to obtain semiclassi- 
cal interpretations of resonance phe nomena. We have 
recently used it in Refs. [15J and [17] in order to under- 
stand, from a new point of view, some aspects of the 
resonant Schwarzschild BH. Our results have permitted 
us to established more particularly, on a rigorous basis 



Appendix A: From Regge poles to complex 
quasinormal frequencies in arbitrary dimensions 

In this appendix, we begin by extracting the resonant 
part of the form factor associated with the S'-matrix de- 
fined in Sec. II. We then interpret it in terms of "surface 
waves", each one associated with a Regge pole. This 
last result permits us to derive the semiclassical formu- 
las (fT9")) - (|21[) which are crucial to construct, in Sec. Ill, 
from the Regge trajectories, the spectrum of the complex 
frequencies corresponding to the weakly damped QNMs. 

Let us consider the scattering of a monochromatic 
scalar plane wave of frequency ui by a d-dimensional static 
and spherically symmetric BH. Without loss of general- 
ity, the corresponding form factor can be written in the 
form 



d-4 

2- _ T fd-3 



2>/7r(zcj) 



x^2{2e + d-3)[S i {u))-l]C i 2 (cos6»i) 



(Al) 

where I is the ordinary angular momentum index, 

2 (z) are the Gegenbauer polynomials (see, for exam- 
ple, Ref. [13|) and Sg(uj) are the diagonal elements of the 
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A Im X 



49 1 



v r : : * 

(d-3)/2-lV jd-3)/2 (d^)/2 + 1 (d-3)/2+2" 



ReX 



We can then deform the path of integration in Eq. (|A2[) 
taking into account the possible singularities (see Fig. [T|) , 
i.e., the poles of the S'-matrix lying in the first quadrant 
of the CAM plane, or in other words, the Regge poles 
A„(oj) with n = 1,2,... By Cauchy's theorem we can 
then extract from Eq. (|A2[) a residue series over Regge 
poles which permits us to obtain the resonant contribu- 
tion in the form 



f F (u,9 1 ) = -2n 



2^ _ r f,l-:\ 



E- 



2-Jlr(iuj) 2 V 2 
A„(w)r n (oj) 



c: 



^sin[7r(A„(a;)-(d-3)/2)] 



Here we have defined 



i_ 3 (-COS( 



(A4) 



FIG. 1: The Watson integration contour. 



S'-matrix defined by Eq. (1171) . By means of a Sommcrfcld- 
Watson transformation [4ll. |42|| (see for example Ref. [2(| 
for a more recent presentation), we can extract in two 
steps the resonant part of the form factor (|A1|) . We first 
replace the discrete sum over the ordinary angular mo- 
mentum I by a contour integral in the complex A-planc. 
We obtain 



i(w) = residue [S 



x-(d-3)/ 2 (w)] A=An(w) . (A5) 



It should be noted that f(u>,9i) differs from /p(w,#i) 
by a background integral over the contour T (see Fig. [1]) 
which does not play any role in resonance phenomena. 

By using the asymptotic expansion (see Ref. [40 ]) 



f(w t 1 )=i 



2^ _ r A/-:$ 



A 



S A _ d^_3 (cj) — 1 



sin[7r(A-(d-3)/2)] 



cos^i) dA. 

(A2) 



Here C is the integration contour in the complex A-plane 
displayed in Fig. [TJ Furthermore, SA-fd-syaO^) is now 
the analytic extension of Se(to) into the complex A-planc 
(CAM plane) which is regular in the vicinity of the posi- 

d-3 

tive real A axis, and C 2 d _ 3 (z) is the analytical extension 
of the Gegenbauer polynomials which is defined by 



c: 



,(*) 



r(A + ^ 



r(d-3)r(A- ^) 



, d — 3 , d— 3d— 2 1 — z\ 

xF[ -A+— ^,A + ^— ) .(A3) 



d-3 



F(A- 



2 J 



e iA(7r-6»i)-i(d-3)7r/4 , e -iA(7r-0i) + i(rf-3)7r/4 



2[tt(A - (d - 3)/2)] 1 /2[( sin 1 )/ 2 ](rf-3)/2 



(A6) 



as |A| — > 00, which is valid for |A|sin6*i > 1, as well as 
the relation 



sm7rz 



+00 

= -2i e l < 2m+1 *> 
m=0 



(A7) 



which is true if Im z > 0, we can write 
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2in 



+ 00 



(2zw) 



— V 

d-2 / y 



r[X»(w) + (d-3)/2] 



n=l 



r [A„(w) + 1/2] [tt(A„H - (d- 3)/2)]V2[( sin ^)/2](rf-3)/2 



^ 1 e iA„(w)(27r-ei+2m7r)-'i(d-3)m7r-i(d-3)7r/4 _|_ e iA„(u))(ei+2m7r)-i(d-3)m7r+i(d-3)7r/4 



(A8) 



In Eq. (jA8|) . exponential terms correspond to "sur- 
face wave" /diffractive contributions. Keeping in mind 
the time dependence exp(— iu>t) and recalling that a 
given Regge pole A„(w) lies in the first quadrant of the 
CAM plane, exp[zA n (u;)(0i)] (resp. exp[?A„(w)(27r — 61)]) 
corresponds to the so-called nth "surface wave" prop- 
agating counterclockwise (resp. clockwise) around the 
BH and Re A„(w) represents its azimuthal propaga- 
tion constant while Im A n (o;) is its damping constant. 
The corresponding exponential decay exp[— Im A„(cj)#i] 
(resp. exp[— Im A„(o;)(27r — 0±)]) is due to continual rera- 
diation of energy. Moreover, in Eq. (|A8I) . the sum over m 
takes into account the multiple circumnavigations of the 
"surface waves" around the BH as well as the associated 
radiation damping. Finally, it should be noted that, in 
Eq. (| A8|> . the presence of the factor exp[— i(d — 3)to7t] 
accounts for the phase advance due to caustics. 

The resonant behavior of the BH can now be un- 
derstood in terms of "surface waves" . Let us consider 
the nth "surface wave", i.e., the "surface wave" de- 
scribed by the Regge pole A„(w). When the quantity 
Re A n (u;) — (d— 3)/2 coincides with an integer, a resonance 
occurs: it is produced by a constructive interference be- 
tween the different components of the surface wave, each 
component corresponding to a different number of cir- 
cumnavigations of the BH [see Eq. (IA8|) ]. the resonance 

excitation frequencies uj^ appearing in the Breit-Wigner 
formula (|18p and generated by the nth "surface wave" 
are therefore obtained from the Bohr-Sommerfeld-type 
quantization condition 



Re A, 



£■ 



d-3 



6 N. 



(A9) 



Now, by assuming that uj is in the neighborhood of cj£ , 

we can expand X n (oS) in a Taylor series about oj^} and 
write 



A„(w) w £ + 



d — 3 dReA n (w) 



+i!mA„ {^fT) + i 



dim A„(w) 



du) 



(A10) 



Then, by replacing (|A10[) in the term 
sin [n(X n (uj) — (d — 3)/2)] of (|A4|) . we can see that 
fp(uj,9i) presents a resonant behavior given by the 



Breit-Wigner formula (TT8"|) with 

Tgn ImA„(oj)[d/cfcjReA n (cj)] 



2 [d/dujRe\ n (uj)} 2 + [d/dujlmX n (uj)} 2 

(All) 

Furthermore, it should be noted that in the fre- 
quency range where the condition \d/du) Re A n (oj)| S> 
\d/duilmX n (u})\ is satisfied, (IA11|) reduces to 



2 



Im Xn(oj) 



d/duj Re X n (ijj) 



jo) 



(A12) 



From the Regge trajectories, i.e., the curves traced out 
in the CAM plane by the functions X n (uj) for lu € [0, +00 [, 
formulas (|A9|) and (|A11I) or (|A12[) permit us to construct 
semiclassically the resonance spectrum of the BH or, in 
other words, the spectrum of the quasinormal complex 
frequencies. Equation (|A9I) provides the location of the 



excitation frequencies ui^ n of these complex resonances 
while Eq. (|A11[) provides the corresponding widths. Of 
course, the reasoning leading to (IA9|) and (|A1 1|) from 
(|A1[) is based on many assumptions that are not nec- 
essarily satisfied in practice. But, in general, it permits 
one to describe rather correctly quasinormal complex fre- 
quencies lying near the real w-axis. 



Appendix B: Regge poles and QNMs of the 
Schwarzschild-de Sitter black hole 

For the d-dimensional Schwarzschild-de Sitter BH 34 
the function f(r) reads 



(Bl) 



Here, Tq is a parameter associated with the mass M of 
the BH by Eq. (|43|) and L is a characteristic length linked 
to the cosmological constant A by the relation 



A = 



(d-l)(d-2) 
2L2 ' 



If 



< 







i) 










L 2 



2 1 (d-3)/(d-l) 



(B2) 

< 1 

(B3) 
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the equation f{r) — has only two real roots and 
rco with < < rco (see, for example, Appendix A of 
Ref. (37j )• They correspond to the radii of the BH horizon 
and of the cosmological horizon. For r E]rh, J~co[, we have 
f(r) > 0. In that case, there also exists a photon sphere 
located at 



and we have 



r c = r 



d-1 



l/(d-3) 



(B4a) 



with r c €]rh,,rco[ and the associated rj c parameter is 
given by 



Vc 



y/d-3. 



(B4b) 



r c and r\ c are both independent of the cosmological con- 
stant and are those already obtained for the ordinary 
Schwarzschild-Tangherlini BH [compare with Eqs. (|44a[) 
and (|44bl) ]. 

Here, it is important to note that, at first sight, the 
formalism we have previously developed does not apply 
to this gravitational background which is not asymptoti- 
cally flat and for which the function f(r) is positive only 
in the finite interval ]r/j, rc [- In particular, the proof of 
the semiclassical formulas (|T9")) - ([2"Tj) given in Appendix A 
is not valid because it is based on the notions of S-matrix 
and form factor which, here, cannot be used. But we can 
easily circumvent these difficulties: 

• Indeed, for f(r) given by (|B1[) the potential Ve(r) 
obtained from ([15)) satisfies the boundary condi- 



tions lim r _>. rh Vg{r) — and lim r 



V t (r) = 0. 



As the tortoise coordinate r* (r) provides a bijection 
from ]r/ l; rc [ to ] — oo, +oo[, it is formally possible 
to define a kind of 5-matrix from the relations (|16[) 
and ([Tf|. 

• The structure of this 5-matrix in the complex em- 
plane and in the complex A-plane allows us to con- 
sider the spectra of the quasinormal frequencies and 
of the Regge poles. 

• It is then possible to establish a semiclassical con- 
nection based on (fT^I) - (f^T]) between these two spec- 
tra. This can be achieved by constructing, from the 
Regge modes, the diffractive part of the Feynman 
propagator associated with the scalar field, i.e., by 
extending a formalism introduced a long time ago 
by Sommerfeld [42| as an alternative to the usual 
approach of scattering [4~fl | considered in Appendix 
A. Such an approach has been used in Ref. [16| for 
the BTZ BH. 



1. The four-dimensional Schwarzschild-de Sitter 
black hole 



For d = 4, the condition (|B3[) implies 

n 2 ™ 2 1 

0< — <1 



(B5) 



r c = 3M, 
Vc = 1, 



(B6a) 
(B6b) 



as well as 



29 137M 2 5 / 27M 2 \ 

' 18 y T?~ ^ a ( n ) 



K = 



216 8L 2 
1357 1505M 2 24765M 4 



15552 288L 2 64L 4 
305 35M 2 705M 4 



3888 72L 2 16i 4 
313 3971M 2 24765M 4 



7776 144L 2 32L 4 
" 235 235M 2 705M 4 

+ 



1944 36L 2 8i 4 

bn 



(B7a) 



a(n) 2 , (B7b) 



a(n) 2 , (B7c) 



27M 2 /(1 - 27M 2 /L 2 ) uj 2 + a n + a{n) 2 ' 

(B7d) 



The WKB approximation (|3T|l for the Regge poles leads 
to 



A„(w) 



27Af 2 



1-27M 2 /L 2 
+ia(n) [1 + e„(cj)] 



w 2 + a n + 2a(n) 2 e n (uj) 



1/2 



(B8) 



and their high frequency behavior (1351) provides 



X n (uj) = 
+ia(n) 



3n/3M 



y/1 - 27 M 2 /L 2 a n 



yJ\~27M 2 jL 2 
(1 - 27M 2 /L 2 )b r , 



27M 2 uj 2 



O 

uj— >-i-oo y oj c 



(B9) 



Here, it is worth noting that, while the radius of the 
photon sphere and the r\ c parameter characterizing the 
instability of the circular null geodesies do not depend 
on the cosmological constant [see formulas (|B6[) ] . this is 
not the case as regards the properties (dispersion rela- 
tion and damping) of the "surface waves" lying close to 
the photon sphere [see Eqs. (|B8[) and (IB9|) ]. As a conse- 
quence, the cosmological constant affects strong gravita- 
tional lensing by the four-dimensional Schwarzschild-de 
Sitter BH. This result corroborates, in a semiclassical 
framework, the recent analysis of Rindler and Ishak [43[ 
concerning the role that the cosmological constant plays 
in the bending of light around a concentrated mass (for 
references on this hot subject, we refer to the bibliogra- 
phy of the recent paper by Ishak, Rindler and Dossett 

The resonance excitation frequencies and the damping 
of the QNMs given by the general formulas (I4U1) reduce 
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to 



complex frequencies. We have 



U) 



(0) 
tn 



Tin 
2 



N /l-27Af 2 /i 2 
3n/3M 

l\ _ On 

2 J 2£ 
V /1-27M 2 /L 2 
3%/3M 



a{n) 



O 



2P 



O 



(BlOa) 
1 

(BlOb) 



and 



A„(w) 



d-3 



1 - 



(d-l)rg 
(d-3)L 2 



Wd- 3a(n) + O (- 



(Bll) 



The leading-order terms of (|B10|) have been obtained 
in Ref. [14J. In Ref. [12j, Dolan and Ottewill have ob- 
tained for n — the expansions of o4? and Yg. n /2 up 
to order l/(£ + 1/2) 4 . Our results (jBlOl) are consistent 
with their Eq. (72). Finally, it should be noted that, 
for L — > +oo, all the results obtained for the Regge 
poles and the complex quasinormal frequencies of the 
four-dimensional Schwarzschild-de Sitter BH reduce to 
the Schwarzschild BH results of Sec. IV.A.l. 



2. Leading-order terms for the Schwarzschild-de 
Sitter black hole of arbitrary dimension 

In the d-dimensional case, we can easily derive the 
leading-order terms of the Regge poles and quasinormal 




d-3 L {d-l)rl ( 1 



- o , 



(B12b) 



Formulas (IB12 | are in agreement with the results ob- 
tained in Ref. [lflj. Furthermore, it should be noted that 
for L — > +oo, from (|B11[) and (|B12[) . we recover the 
Schwarzschild- Tangherlini BH results of Sec. IV. A. 4. 
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